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It has previously been established that the logarithmic-depth approximate quan- 
tum Fourier transform (AQFT) provides a suitable replacement for the regular QFT 
in many quantum algorithms. Since the AQFT is less accurate by definition, poly- 
nomially many more applications of the AQFT are required to achieve the original 
accuracy. However, in many quantum algorithms, the smaller size of the AQFT 
circuit yields a net improvement over using the QFT. 
t*"*- ■ This paper presents a more thorough analysis of the AQFT circuit, which leads 

to the surprising conclusion that for sufficiently large input sizes, the difference 
between the QFT and the logarithmic-depth AQFT is negligible. In effect, the 
AQFT can be used as an direct replacement for the QFT, yielding improvements 
in any application which does not require exact quantum computation. 

& . 

The quantum Fourier transform is a fundamental component of many of the significant 
quantum algorithms known to date. The main concern with the efficiency of the QFT 
algorithm is its use of controlled phase-shift gates which involve very small phases. As 
the number of qubits in the input, n, increases, the QFT algorithm requires exponentially 
smaller phase shifts, which may be increasingly difficult or even infeasible to implement 
physically. 

Fortunately, for applications using the QFT which do not require an exact answer, the 
idea of an approximate QFT was suggested[3], in which the phase shift gates requiring 
the most precision are omitted, since these are also the phase shift gates which have the 
least effect on the output. There is clearly a trade-off between the accuracy of the AQFT 
and its circuit complexity. Barenco, Ekert, Suominen and Torma established a lower 
bound for the probability that the AQFT returns the correct answer, and concluded that 
the accuracy of the QFT can be achieved by 0(n 3 /m 3 ) iterations of the AQFT, where 
n is the size of the input register, and m is the depth of the AQFT circuit. Usually, we 
have m = O(logn). 



1 Introduction 
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Figure 1: A quantum circuit that performs a quantum Fourier transform 

However, using a more thorough analysis of the AQFT circuit, we derive a much 
better lower bound on the accuracy of the AQFT which allows us to conclude that the 
accuracy of the QFT can be achieved by just 1 + o(l) iterations of the AQFT. 



2 The Quantum Fourier Transform 

The quantum Fourier transform is a quantum circuit which performs a discrete Fourier 
transform on the complex-valued vector of 2 n probability amplitudes associated with an 
n-qubit quantum system. Specifically, given an n-qubit state as a superposition of basis 
states |0), |1), . . ., |2 n — 1), the QFT maps each basis state \j) to 
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Figure 1 gives a quantum circuit which performs a QFT. The input register contains 
an n-qubit basis state \x) = \x\X2 - ■ ■ x n ) . The gates labelled H are Hadamard gates, 
and the gates labelled R m represent a series of one-qubit phase rotation gates. For each 
integer m > 2, the gate R m shifts the phase of the |1) component of the input qubit by 
a factor of e 2 ""*/ 2 " 1 . Note that the output state of the QFT can be factored into a tensor 
product as indicated in the diagram. 

In practice, the QFT is usually used in reverse, where the states 

_L ( | 0)+e 2^-V)| l)) 

are prepared for p = 1,2, ... ,n by some other means, and used as input into the in- 
verse QFT circuit in order to obtain an estimate of the value of from the state 
\x) = \x1X2 ■ ■ -x n ). Specifically, is estimated by = (O.X1X2 . . .x n ), where the paren- 
theses denote a binary (base-2) fraction. The estimate is the nearest multiple of 2~ n 
to with probability and one of the two nearest multiples with probability 

As for the generation of the required input states, since it is not important how they 
are generated, we may assume that they are generated by using the state ^=(|0) + |1)) 

as the control qubit for a controlled- U 2P 1 gate, where U is a quantum gate with as an 
eigenvalue. 
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Figure 2: An individual trial 



Griffiths and Niu observed that it is not necessary for the entire inverse QFT to be 
performed at once [4]. Instead, since each output qubit, \x p ) is controlled only by higher- 
indexed qubits, and the last qubit \x n ) can be obtained independently of all other qubits, 
they suggested a "semiclassical" version of the inverse QFT circuit in which each output 
qubit is obtained one at time, starting with \x n ). Proceeding one qubit at a time, the 
results from measuring the previous qubits are used to classically determine whether a 
particular phase rotation gate will be applied to the current qubit. This process does not 
affect the final answer, but has the advantage of replacing the controlled-rotation gates 
in the inverse QFT with regular rotation gates. Thus, we can analyze the entire inverse 
QFT circuit simply by looking at each individual "trial". 

In the circuit in Figure 2, the controlled rotation gates are replaced with a rotation 
by a phase of — Xp, which is determined by the results of the previous trials, from \x n ) to 
\x p+ i). Specifically, Xp — (0.0x p+ i . . . x n ). 



3 Approximate QFT Circuits 

In Coppersmith's AQFT circuits[3], instead of performing all phase rotations R p that 
are used in the QFT algorithm, we set a lower limit on the amount of phase shifted by 
any particular gate and ignore any phase rotation gates that do not shift by at least 
that amount. In other words, given a positive integer threshold, m, we will define the 
approximate QFT circuit AQFT m to be the circuit formed from the QFT, except that 
phase rotation gates R p are ignored whenever we have m > p. Note that AQFT n is simply 
the regular QFT circuit. As before, the inverse AQFT can be divided into individual 
trials. 

Suppose we are given a phase 0, and we are to use the AQFT m circuit to find an 
approximation = (O.X1X2 ■ ■ -x n ). The individual trial for bit x p starts with the input 
state |0) + e 2m( - 2P w |1), and performs a phase rotation of e~ 2mXp on the probability 
amplitude of |1). Here, Xp corresponds to the phase rotation amount performed by the 
regular QFT circuit, (0.0£j + i£j + 2 . . -x n ), but without the effect of any gates which are 
smaller than the given threshold m. Therefore, we have Xp = (0.0x p+ i . . . x p+m _i) when 
p < n — m, and Xp — (0.0x p+ i . . . x n ) otherwise. 

We will now compute the probability that is the nearest integer multiple of l/2 n 
to 0, i.e., |0 — 0| < 2"™" 1 . First, using the double angle formula and induction, we can 
easily establish the fact that for any integer n > 1, we have 

n 2 9 ( sin6> \ 
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Figure 3: A quantum circuit implementing the AQFT 

Now observe that, in general, applying a Hadamard gate to the state |0) -\- e 2nl ((°- x p)+ 5 p) 
and measuring yields \x p ) with probability cos 2 {tt5 p ). However, for the individual AQFT m 
trial for bit x p , we know what S p has to be. We have two cases. 

In the first case, where p + m — 1 < n and \p = (0-0x p+ i . . . x p+m _i), 5 P contains the 
phase rotation amount that was not applied by the AQFT m trial that would have been 
applied by the corresponding QFT trial, and the difference between <p and its nearest 
integer estimate, amplified by a factor of 2 P by the AQFT, n circuit. In other words, we 
have 



where 5 



5 p 2 (O.Xp-|_ m Xp+ m _|_i . . . x n ) -|- IP S, 

Also, if is the nearest estimate of 0, we have 
_ 2 p-n-i < 2 p~i s < 2 p-n-l_ 



Now, since (0.x p+m . . . x n ) reaches its minimum value when all the bits are 0, and its 
maximum value when all its bits are 1, we have 
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Adding the two inequalities together, we obtain — 2~ m < —1~ n ~ < 5 P < 2~ m . This 
means that final measurement will yield \x p ) with a probability of at least cos 2 (7r2 _m ) in 
this case. 

In the second case, where p > n — m+1, we simply have 5 P = 2 P " 1 5. Since \5\ < 2~ n ~ 1 , 
we have an error probability lower bound of cos 2 (|2 p_n ~ 1 ). The probability P of getting 
the correct output with the AQFT m circuit is simply the product of each individual bit 
trial being correct. Putting these lower bounds together, we can give a lower bound for 
P: 
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It is clear that if m grows too slowly in comparison to n, then the expression will 
approach asymptotically as n approaches oo, making this bound useless. If we take 
Tt > log 2 n + 2 so that 2 m > 4n, we would have 

Hwr- > ± p (£))", 

giving us 
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Since Cleve and Watrous[2] establish a lower bound of f2(log n) on the depth of any AQFT 
circuit with a constant error bound, requiring m > log 2 n + 2 is a reasonable restriction 
on a logarithmic-depth AQFT circuit. 

Asymptotically, for m > log 2 n + 2, the lower bound for the accuracy of the AQFT m 
circuit approaches the lower bound for the accuracy of the full QFT, meaning that for 
sufficiently large n, the QFT may be substituted by the AQFT m circuit with only a 
negligible effect. This is a significant improvement, since the full QFT on an n-qubit 
register requires 0(n 2 ) gates to implement, while the AQFT m circuit uses only 0(nm) = 
O (n log 2 n) gates to achieve nearly the same level of accuracy. 

Finally, we will also establish a lower bound for P when we have m > log 2 n + 2 for 
some fixed value of n, since an asymptotic bound gives no information about the value 
of P for practical values of n. We have 
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using the Bernoulli inequality, which states that (1 + t) c > 1 + ct whenever t > — 1 and 
c > 0. Since x > sinx when x > 0, we have 
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This indicates that even for smaller values of n, the probability of the AQFT m circuit 
returning the best estimate of 4> is high enough to make it a practical alternative to using 
the full QFT circuit. 

Note that since we stipulated that m > log 2 n + 2, this result is only meaningful if 
n > A. Otherwise, the logarithmic-depth AQFT m circuit is simply the full QFT, with 
lower bound P > \. We can now establish the fixed bound 
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for AQFT m circuits where m > log 2 n + 2. 
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4 Conclusion 



The fixed bound of P > \ — derived here is much stronger than the bound of 



previously established by Barenco et al. in [1]. They conclude that 0(n 3 /m 3 ) iterations 
are necessary which means that the AQFT m yields an advantage only when the QFT 
itself is affected in the presence of decoherence, as the AQFT. m would be less susceptible. 

The new analysis presented in this paper establishes the effectiveness of the AQFT m 
as a direct replacement for the QFT for a suitable size of input register. Furthermore, 
since we have 



when m is fixed at log 2 n + 2, it was not previously known that the AQFT m circuit would 
perform so well for large n. This is a significant result for the feasibility of any large-scale 
quantum computation requiring the use of the QFT circuit. 
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